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ABSTRACT 

Molecular hydrogen (H2) absorption features observed in the line-of-sight to 
Q2348— Oil at redshift Zabs —2.426 are analysed for the purpose of detecting a pos- 
sible variation of the proton-to-electron mass ratio fi = mp/m^. By its structure 
Q2348— Oil is the most complex analysed H2 absorption system at high redshift so 
far, featuring at least seven distinctly visible molecular velocity components. The mul- 
tiple velocity components associated with each transition of H2 were modeled simul- 
taneously by means of a comprehensive fitting method. The fiducial model resulted in 
A/z//x = (— 0.68± 2.78) x 10~^, showing no sign that ^ in this particular absorber is 
different from its current laboratory value. Although not as tight a constraint as other 
absorbers have recently provided, this result is consistent with the results from all 
previously analysed H2-bearing sight-lines. Combining all such measurements yields a 
constraint of |A/i//i| ^ 10^^ for the redshift range z = (2-3). 

Key words: methods: data analysis, quasars: absorption lines, cosmology: observa- 
tions 



1 INTRODUCTION 

The spectrum of molecular hydrogen is known to be a test- 
ing ground to search for temporal variation of the proton- 
electron mass ratio, fi = mp/m^, on cosmological time- 
scales (Thompson 19751. The spectral lines of the Lyman 



and Werner absorption bands of H2 shift into an atmo- 
spheric transmission window (A > 3050 A) for absorption 
systems at redshift > 2, and thus become observable with 
ground-based optical telescopes. Wavelengths of H2 spec- 
tral lines as observed in high redshift absorption systems 
are compared to the wavelengths measured in the labora- 



tory for which now very accurate calibrations exist (Philip 
et al.|2004[[Salumbides ct al. 2008 ,'B ailly et al.|2010|) as well 



as for the HD isotopologue (Ivanov et al.|2008| [Ivanov et al.| 
2010), which is currently also detected at high redshift. Val- 
ues for the sensitivity coefficients, Ki, expressing the shift 
of each spectral line as a result of a drifting /x, have been es- 
tablished to sufficient accuracy for H2 ( [Meshkov et al.|2006| 



Ubachs et al.|2007 1 . Hence a condition is accomplished that 
searches for cosmological fj, variation via H2 spectra solely 
depend on the accuracy of the astrophysical data. 

From the large number of damped Lyman-Q systems 
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(DLAs) identified only some 20 are known to harbor de- 
tectable H2 spectral featu res ([Ledoux e t al. '2003| [Srianand] 
et al.||2005| [Ubachs et al.||2007| ). As" was discussed recently 



(Ubachs et al. 20111 only a handful of those high redshift 
absorption systems can deliver an H2 spectrum of suffi- 
cient quality, i.e. a reasonable signal-to-noise (S/N) ratio 
on the continuum level to be obtained in reasonable data 
collection times at large-dish telescopes, and a large num- 
ber of H2 transitions (> 50). The high quality spectra of the 
Q0347-383 and Q0405-443 systems obtained with the ESQ 
Very Large Telescope (VLT), equipped with the high reso- 
lution Ultraviolet and Visible Echelle Spectrograph (UVES) 
initially yielded an indication for a possible drift in ji, which 
is expressed via A^/jj. = (fiz — ^tiab)/Miab (Reinhold et al 



20061. A more sophisticated re-analysis of the same spec- 
tra by the so-called comprehensive fitting method, also in- 
voking an improved thorium-argon (Th-Ar) calibration, re- 
duced the initial finding of a 4-(t effect of a positive Afi/ ^ 
to a 1.5-2.0-cr effect ( [King et al.][2008[). The same spectra 
were also re-analysed by Wendt & Reimers ( 2008 1 and by 



Thompson et al. (20091, while Wendt & Molaro (20101 re- 



observed Q0347-383 with VLT-UVES to find a 1.5-cr effect. 
The highest quality system observed so far is J2123— 0050. 



Observations with both the Keck Telescope (Malec et al. 
•2010) and VLT ([van Weerdenburg et al.][2011[) were incor- 



porated in a /i-variation analysis. The two studies yielded 
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constraints on Ajj,/ fi which are in very good agreement 
with each other. The averaged constraint for J2123— 0050 
is A/i//x = (7.6 ± 3.5) X 10"^ The system Q0528-250 has 
been at the heart of p-variation analyses since the first at- 



using the object profile weights from its corresponding 
quasar exposure. Further details of the wavelength calibra- 



tempts by Varshalovich & Levshakov ( 1993 \ to derive a con- 



straint on AfJ,/ fJ, from a low-resolution spectrum. King et al 



( 2008 1 derived a tight constraint from a comprehensive fit- 



ting analysis of a VLT-spectrum of Q0528— 250; a recent 
re-observation of this high-quality system, again at VLT, 
confirmed the previous conclusion of a tight constraint on 
A/i//x ( King et al.|20H I. Based on the combined results de- 
rived from high redshift H2 absorbers it can be concluded 
that \Afj,/fi\ < 1 X lO"'^ for redshifts z = 2 - 3. 

The present study presents a detailed analysis of the 
search for /i-variation of an absorber system toward the 
quasar Q2348— Oil. The physical conditions of this system 
had been investigated before based on data from the VLT 



(Petitjean et al. 2006 Ledoux et al. 2006 Noterdaeme et 
|al.|2007| ). Q2"348-011 was shown to be an exceptional sight- 
line, with a number of DLAs and sub-DLAs present, and a 
very complex velocity structure in the major H2 absorbing 
system at z ~ 2.426. At least seven distinctly visible absorp- 
tion features are associated with each H2 transition in this 
system, while in other quasar sight-lines H2 absorption (if 
detected) exhibits much simpler profiles consisting typically 
of one or two features. For the present study dedicated ob- 
servations were performed at VLT focusing on improvement 
of S/N and on wavelength calibration of the spectrum. A 
preliminary /.i-variation analysis of this absorber was pre- 
sented in Ubachs et al. (2011). The current work improves 
upon the analysis but the constraint on A/fi is of similar 
precision as presented previously. It is demonstrated that a 
complex velocity structure consisting of at least seven dis- 
tinct H2 velocity features can be disentangled by the method 
of comprehensive fitting ( King et al.|2008 Malec et al.|20ld 
King et al.|2"0TT| ). 



2 DATA 

The new spectrum of the quasar Q2348— Oil was obtained 
with the UVES spectrograph on the ESO VLT on four con- 
secutive nights (2008 August 18-21). Q2348-011 is not a 
particularly bright object (R = 18.31) so a long exposure 
time was required. Fifteen exposures were taken, making 
up a combined total of 19.25 hours of observation. After 
each science exposure, without any intervening grating re- 
sets, a Th-Ar exposure was recorded for calibration pur- 
poses. Seeing was in the range between 0'.'53 and 2'.'32. Be- 
sides the exposures noted above, the final spectrum incorpo- 
rates data obtained from the ESO data archive (observations 
from 2003, October 29-30). They contribute additional 4.50 
hours but do not have the individually taken Th-Ar calibra- 
tion frames. 

The raw 2D exposures were bias-corrected, flat-fielded, 
and the quasar flux extracted using the Common Pipeline 
Language version of the UVES pipeline^ The wavelength 
calibration was established by extracting each Th-Ar frame 



tion procedure are described in Murphy et al. (20071. The 
root-mean-square wavelength calibration residuals were typ- 
ically ~70 m s^^ in the Lyman-a forest portion of the 
spectrum, where the H2 lines fall. The wavelength scale of 
each quasar exposure was corrected from air to vacuum, 
and to the heliocentric reference frame, and the flux was 
re-dispersed on to a common log-linear scale before being 
co-added, using UVES_POPLEIij^ a code speciflcally written 
to re-disperse and combine reduced exposures from UVES. 

The spectrum of Q2348— Oil covers wavelengths 3572- 
9467 A, with gaps at 4520-4621 and 7505-7665 A. The 
final spectrum has a resolving power of R ~ 57,000 and 
63,000, and the log-linear dispersion is set to 2.5 km s~^ and 
1.5 km s"^ in the blue (< 4000 A) and in the red (> 4000 A) 
parts, respectively. CCD pixels were binned by a factor of 2 
in both spatial and spectral directions for all exposures in 
the blue, and no binning was applied for the red exposures. 
In the blue part where H2 features are detected the average 
S/N ratio is about 25. 



3 THE COMPREHENSIVE FITTING METHOD 

The main aim of the present analysis is to determine the 
spectral positions of the H2 absorption features present in a 
quasar spectrum as accurately as possible since their overall 
pattern defines the value of A/i//i: each absorption line of the 
H2 Lyman and Werner bands shows a unique shift for a given 
change of fj, that depends on the vibrational and rotational 
quantum numbers of the upper and lower energy states. The 
shifted wavelength Xi is related to the rest wavelength Ao by 

X,^Xo{l + z){l + K,^) (1) 

where z is the redshift of the absorber, and Ki is the sensi- 
tivity coefficient, different for each line. To achieve this goal 
a so-called comprehensive fitting method can be applied. 
This method requires fitting the absorption profiles of mul- 
tiple H2 transitions simultaneously. It relies on a physical 
assumption that all transitions arise from the same cloud 
(collection of clouds) of molecular hydrogen and, therefore, 
they share the parameters describing the properties of the 
clouds (see e.g., [Malec et al.|2010[ ). 

Each H2 transition can be detected in one or several 



absorption features depending on how many molecular hy- 
drogen clouds are penetrated by the light of the quasar. We 
refer to each fitted component as a velocity component, as 
they are situated closely in velocity space. It is assumed that 
a given velocity component has the same redshift, ^abs, and 
Doppler width, 6, in all transitions independent of the ro- 
tational level, J. Each J level has a different population, 
represented by the column density A'^, but for a given ve- 
locity component it is the same for all transitions from the 
same J level. To analyse the spectrum we use the Voigt pro- 
file fitting program VPFn]^ which permits to tie the free 



^ |http://w ww.eso.org/observing/dfo/quality/UVES/pipeIine/pipe| 
_reduc.html. 



^ Developed and maintained by M. Murphy; 
[http:/ / astronomy.swin.edu.au/ ^mmurphy/UVES-popler 
EFtp : / / www. ast . cam .ac.uk/~rfc/vpiit. html. 
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Table 1. A list of all the confirmed detections of metal absorbers 
in the spectrum of Q2348— Oil. Species whose transitions overlap 
some of the fitted H2 regions are in bold. 



Redshift 


Species 


0.774 


Mg 11 


0.863 


Mg II, Ca II 


1.444 


C IV, Mg II 


2.426 


C IV, Si IV, N V, C I, C I*, P II, 




Zn II, Al III, Fe II, S II, Si II 


2.582 


C IV, Si IV, Si II 


2.615 


C IV, Si IV, Fe II, Ni II, Cr II, Al II, 




Si II, I, N I 


2.739 


C IV, I 


2.929 


C IV, Si IV 


2.969 


C IV, N V 



parameters between various transitions. The optimal model 
is found by fitting all suitable transitions at the same time. 

The differential shifts of the molecular hydrogen lines 
are only known to be possible through the A^/ jj, param- 
eter (see Eq. It is added as a free parameter after all 
other parameters have been optimized. The comprehensive 
fitting approach allows minimization of the number of free 
parameters in the fit and, thus, the reliability of the absorp- 
tion model which eventually translates to the robustness of 
Afi/jj, resulting from the model. 



4 ANALYSIS 

4.1 The absorption system 

The wavelengths covered by the VLT-UVES spectrum pro- 
vide 58 H2 transitions for rotational levels J = — 5. All the 
observed lines arise from the Lyman band. In Fig.[T]a couple 
of H2 transitions are displayed. No lines of HD were detected 
in this system. In the ideal case this spectrum would pro- 
vide a sample of 58 x 7 molecular hydrogen absorption lines. 
However, not all of them are suitable for the analysis. First 
of all, the broad Ly-/3 of an additional DLA at 2: ~ 2.62 
falls in the wavelength range where II2 lines are detected, 
namely near 3710 A. The break at ^3590 A might be pro- 
duced by another absorber at z ~ 2.93. These features make 
some of the H2 transitions unavailable, and strongly damp 
some others. Also, the weaker H I absorbers of the Ly-a 
forest contaminate many of the regions with H2 lines. Gen- 
erally, the sight-line of the Q2348— Oil is rich in absorbing 
systems. The full list of detected metal absorbers is given 
in Table [l] Some of the absorbers have their metal absorp- 
tion lines overlapped with II2 lines. The metal absorption 
profiles are constrained by fitting their counterparts in the 
red part, e.g., Fe II A1096 of the z ~ 2.426 absorber can be 
constrained by tying it with Fe II A 1608. 

Since the velocity structure of H2 consists of at least 
seven velocity features spread over ~ 250 km s~^, nearby 
lying transitions exhibit mutual overlap of velocity compo- 
nents. The self-blending of the H2 transitions was often the 
decisive factor when setting the boundaries of usable fit- 
ting regions (i.e. regions over which was calculated): only 
rarely did they include just one transition (see Fig. [5] as an 
example) . 



B 1 



o 
Z5 




-150 -75 75 
Relative velocity (km s~^) 

Figure 1. Two H2 and two Fe II transitions are displayed on a 
common velocity scale, centered at zq = 2.4263230. In the H2 
velocity profile 7 distinct features are visible, while the Fe II pro- 
files consist of at least 14. Note that the third velocity component 
from the left in the H2 LORl profile, and the right-most compo- 
nent in the H2 L1P2 profile are blended with other intervening 
transitions in Ly-a forest. Due to such blends, the signal of at 
least 7 absorption components is evident only when multiple H2 
transitions are explored simultaneously. 



Twenty-three regions with 32 different molecular hydro- 
gen transitions and eight regions with relevant metal absorp- 
tion profiles were selected for further analysis. Although the 
total number of H2 transitions available in this spectrum is 
reduced due to blends with H I lines, the scarcity of infor- 
mation is at least partially compensated for by the presence 
of multiple velocity components associated with each tran- 
sition. Altogether, each transition contributes a higher in- 
formation content when it is imprinted on several velocity 
components. A part of the spectrum with all available and 
all fitted H2 transitions indicated is provided in Fig. [3] 



4.2 Exploring the velocity structure of H2 

Seven absorption features can be distinguished by eye in the 
velocity structure of molecular hydrogen in the Q2348— Oil 
spectrum (see Fig. [l] a lso [Noterdaeme et al.| ([2007)) . How- 
ever, as it was shown by |Murphy, Webb fc Flambaum ( 2008 1, 
a single spectral feature might need to be modeled by more 
than one velocity component so it is plausible that in the 
case of Q2348— Oil more than seven velocity components are 



4 



J. Bagdonaite et al. 




3765 3766 3767 3768 3769 
Wavelength (A) 



3770 



3771 



Figure 2. H2 absorption profiles of two transitions: L1P3, and 
L1R4, resulting from the fiducial model. Since the absorption 
arises in multiple widely spread components (indicated by the 
red tick marks for L1P3, and by the blue tick marks for L1R4), 
the two profiles overlap. The upper graph shows the normalized 
residuals of the fit; the dashed lines mark ±1 cr limits. 



actually required to produce a model which is both physi- 
cally realistic and statistically acceptable. If a single Voigt 
profile is fitted to a feature which seems single but another 
weaker blending line is actually present, then the centroid 
returned from this fit will be shifted towards the blending 
line. The centroid wavelength of the fitted feature might be 
expressed as the intensity-weighted mean of the two blended 
lines. The positions of the lines play a crucial role in /i- 
variation analysis and this kind of underfitting may affect 
the results and therefore should be avoided. 

The approach to reaching the optimal model of H2 ab- 
sorption is to keep adding velocity components until a statis- 



tically acceptable fit is achieved (e.g., Malec et al. 2010 King 
et al.|20lT l. The initial absorption model included seven ve- 
locity components of H2. Voigt profile parameters {z, b, N) 
of corresponding components are tied as described in Sec- 
tionjs] Each of the included regions was revisited and refined 
many times until an adequate fit was reached. By adequate 
we mean that the normalised residuals in each of the fitted 
regions are not too large or too small (i.e. they are dis- 
tributed around [— l-cril-cr]), and no long range correlations 
are seen in them. Finally, after the fit was optimized, a value 
for Afi/ 11 was determined. 

In a subsequent step an additional component was 
added to the left side of the absorption profile (at ~ —150 
km s~^; see Fig. [l]). The reduced chi-square, Xv, of this 
model is slightly better than that of the initial one with 7 
velocity components. The procedure of adding more compo- 
nents was continued further until they started to be rejected 
from most J levels (>12 component models). We note that 
several components get rejected from J — 4 and 5 transi- 
tions in all models. Generally, it is not urmsual that weak 
components are omitted from the lowly populated J levels. 
The higher J levels are less populated than the lower J lev- 
els, so weak components of J = 4 and 5 transitions may not 
be detected. In addition, the Q2348— Oil spectrum provides 
only a few useful regions with J — 4 and 5 transitions, which 
is not enough to amplify the signal of the weak components. 

It should be emphasized that Afi/ jj, is always intro- 
duced as a single additional parameter after all the other 
free parameters of the given model have been optimized. 
This ensures that Afi/fi is not biased away from zero in 
the process of determining which velocity structure is sta- 



1.25 
1.225 
1.2 



7 8 9 10 11 12 

Number of fitted velocity components 

Figure 4. Lower graph: A/i/n values as determined from various 
models. Error bars indicate estimated l-tr statistical uncertain- 
ties. Upper graph: xi of the fit decreases as more components are 
added. A model with 12 components has the lowest xi- Adding 
more components caused inconsistencies in the fit. 



tistically preferred. On the other hand, it could mean that 
Afj./ jj. might be slightly biased towards zero because the pre- 
ferred velocity structure might have 'fitted away' some of the 
relative shifts between transitions which a non-zero Afi/ ^ 
would have caused. However, given that many H2 transi- 
tions are fitted simultaneously, the bias towards zero will be 
very weak. 

Results of the various fitting runs are displayed in Fig.[4j 
together with the goodness-of-fit measure, xt- As it can be 
seen in Fig. [4] the xt monotonically decreases when more 
components are added. For models with 10, 11 and 12 com- 
ponents the difference in Xv is very small though. This is 
also true for the derived Afi/n values - they match well 
within 1-cr. For further analysis the model with 12 velocity 
components has been adopted as fiducial. 



4.3 Robustness of the model 

When all available transitions are included in the 12 compo- 
nent model, the fit delivers Afj./fj. = (-0.68±2.70stat) x 10"^. 
Next the fiducial model is tested for its robustness. In the 
process of creating the model, many assumptions and deci- 
sions have been made, e.g., how parameters are tied between 
transitions and which regions are fitted. It can be tested how 
sensitive the value of A/i//i determined from the model is 
to certain choices made. A description of several performed 
tests is given below, while the outcome is summarized in 
Tabled 

(i) The extended structure of the H2 absorption profile 
can be divided into two parts. As it can be seen from 
velocity plots, the dividing line naturally falls at around 
km s~^ (Fig. [l]). When only the left part of the profile, 
which includes 4 H2 velocity components, was used in the 
fitting, it resulted in A/x/^ = (3.55 ± 4.34stat) x 10"^ The 
right part with 8 velocity components delivered Afi/fj, — 
(-3.87 ± 3.43stat) X 10"^. The larger uncertainty of the left 
side result can be explained by the fact that several of the 
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Figure 3. The part of the Q2348— Oil spectrum containing the H2 transitions. Transitions from different J levels are labeled in different 
colours. The associated ticks show positions of the 7 (visually distinguishable) velocity components. The broad absorption feature at 
3710 A is the Lyman-/3 line of the H I absorber at z ~ 2.62. The break at the shortest wavelengths is due to another H I absorber at 
z ~ 2.93. The green line shows the fitted regions. 
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Table 2. Resulting Afi/fi values and their statistical uncertainties from various consistency tests 
performed on the 12 velocity component absorption model as described in Section |4.3| The b 
F[J] and the b = F[J] models are different in that the former has b parameters tied in all J levels 
of corresponding velocity components whilst in the latter they are only tied within every J level 
of corresponding velocity components. The last line of the table refers to a test, in which only the 
best exposures (calibration-wise) were included in the spectrum. 





b 


Transitions 


H-trans 


Components 


Afi/fi [xlO-^] 






all 


32 


all 


-0.68±2.70 


i 




all 


28 


left 


3.55±4.34 


i 




all 


22 


right 


-3.87±3.43 


ii 




J = [0, 1] 


11 


all 


-0.86±3.50 


ii 




J = [2 - 5] 


21 


all 


-0.46±4.06 


iii 


= F[J] 


all 


32 


all 


-1.20±3.03 


iv 


7^m 


all 


32 


all 


0.68±3.45 



fitted transitions were saturated, thus their centroids are 
determined less accurately. Altogether, the two results are 
consistent within the combined 1.3-a. 

(ii) Instead of allowing all transitions to contribute to a 
single Afj,/^ value, within VP FIT it is possible to fit a differ- 
ent value of Afi/iJ, for each J-level or for J-levels grouped 
in some way. This kind of test allows to quantify the rela- 
tive contributions different J-levels make to the final result. 
As not many transitions are available in the Q2348— Oil ab- 
sorber they were not investigated level by level. |Ubachs "it] 
al. ( 2007 1 have suggested that H2 transitions can be divided 



4.4 Systematic errors 

The statistical error of the fiducial result is larger than ob- 



into a J G [0, 1] set (cold states) and a J 2 set (warm 
states) to examine the impact of temperature: due to the 
para-ortho distribution of H2 the J — 1 state is significantly 
populated even at low temperatures. A test was performed 
where only transitions from the cold states or only those 
from the warm states were used to determine Afj,/fi. The 
two values of Afi/fi match within the uncertainty. It means 
that the two groups of transitions contribute similarly to the 
final combined value of An/fi. 

(iii) All previously described tests were performed on ab- 
sorption models which are based on the assumption that 
corresponding velocity components in all transitions have 
the same 6-parameter, independent of rotational quantum 
state (6 7^ F[J])- This assumption is relaxed to test what 
its effect is on Afi/fi, i.e. in this run different J transitions 
are allowed to have different 6-parameters in corresponding 
velocity components (6 = ^'[J]). No significant discrepancy 
is found between the values determined from the b F[J] 
and the b = F[J] models. 

(iv) The spectrum used in the analysis includes exposures 
from the older 2003 ESO-archive data set which do not have 
the individually taken Th-Ar calibration frames. They make 
up 4.50 hrs compared to 19.25 hrs of more accurately cali- 
brated data. In order to make sure that the final result is not 
affected by the additional (possibly less accurate) data, they 
were excluded in one of the fit optimizations. The newly de- 
termined value is A^i/fi =(0.68±3.45stat) x lO"'^. The result 
is consistent with the one from the primary fit. The uncer- 
tainty is higher since the S/N of the spectrum used in this 
test is lower. 

The performed tests show that the final result on Ap//i 
determined from the 12 component model is not affected by 
a subclass of transitions, velocity components or data. 



tained in other recent studies of H2 studies (e.g., |Malec et al^ 
2010 King et al.||2011 l. In those studies, systematic errors 
were found to be smaller than, but comparable in magnitude 
to, the statistical errors. Thus, while it is still important to 
consider systematic errors in the present study, they are not 
expected to dominate the error budget. 

Every echelle order included in the new Q2348— Oil 
spectrum is calibrated using about 10 Th-Ar lines. The ef- 
fect of the errors of individual Th-Ar lines is small (about 
70 m s~^) and random, so it will average out when many H2 
transitions spread over many orders are used. |Murphy et al.| 
(2007) have shown that possible systematic patterns in the 
Th-Ar calibration residuals result in a deviation of 30 m s~^ 
at most. The Ki values used in the present analysis are in 
the range from —0.015 to 0.018. The effect of the system- 
atic errors can be expressed via: 5{Afi/^) = (Av/c)/ AKi. A 
shift of 30 m s~^ translates to S{Afi/fi) = 0.3 x 10~^ which 
is about 15 per cent of the statistical uncertainty. 

In a recent study targeted at the possible miscalibration 
of the wavelength scale of the UVES spectrograph it is re- 
ported that intra-order distortions of about 200 m s~^ in size 



may be present ( Whitmore, Murphy fc Griest||20To I. Their 
influence can be roughly estimated by dividing the distor- 
tion by the square root of the number of the H2 transitions 
included in the fi t : 200/ \/32 ^ 35 m s~\ 

Mal ec et al.| ( |2010[ ) have considered additional sources 
of possible systematic errors but they are all of the same 
order as the Th-Ar wavelength calibration uncertainties and 
are not expected to have substantial impact on the result of 
the present analysis. It is conservatively assumed that the 
systematic error is at most 25 per cent of the statistical error 
of Afi/n derived in the present analysis. 



5 RESULT 

Based on the results provided in Table [2j 

A^/^i = (-0.68 ± 2.70stat ± 0.66sys) x 10"^ 

is adopted as the fiducial result for the analysis of molec- 
ular hydrogen absorption in the Q2348— Oil spectrum (or 
An/jj. — (—0.68 ± 2.78) x lO^'', adding the uncertainties in 
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quadrature). The result is consistent with no change in fi at 
the 10"^ level. 

The analysis of H2 absorption in the Q2348— Oil spec- 
trum yields a result on A^/jj, which is less tight compared to 
the constraints from previous analyses (Ubachs et al. 2011). 
One of the major causes that led to the larger uncertainty 
is the relatively low S/N ratio of the spectrum. The sec- 
ond cause is that the number of lines is relatively small: 32 
compared to e.g., 90 in the sight-line of J2123— 0050, which 
yields an order of magnitude more precise result on A^//i 
( Malec et al.|2010[[van Weerdenburg et al.|20"Tl| . Part of the 
H2 spectrum in Q2348— Oil is obscured by the neutral hy- 
drogen features of the additional strong DLA at z ~ 2.615 
making some of the relevant transitions unavailable. The 
H2 transitions falling at the short wavelength range of the 
spectrum are especially useful in ^-variation analysis, since 
they exhibit the larger sensitivity coefficients. In the spec- 
trum of Q2348— Oil they are not detected due to the H I 
absorption produced by the additional DLA and another 
strong absorber at 2: ~ 2.93. However, the present analysis 
demonstrates that a complex absorption structure can be 
successfully modelled by means of the comprehensive fitting 
method. 



6 CONCLUSION 

Molecular hydrogen features present at redshift Zabs — 2.426 
in the line-of-sight to the quasar Q2348— Oil were analysed 
to detect a possible variation of the proton-to-electron mass 
ratio on a cosmological time-scale. The constraint derived 
in the analysis is Afi/^i = (-0.68 ± 2.78) x lO"'^, show- 
ing no indication that jj, in this particular absorber is dif- 
ferent from fj, measured in the laboratory. Although being 
less accurate than other recent H2 constraints on A^/ji, 
the result is consistent with those results, which show that 
|A/i//i| < 1 X lO"'^ at z = 2-3. 

In the sample of known high-redshift II2 absorbers, the 
structure of the absorber in the Q2348— Oil spectrum is the 
most complex as it has 7 visually distinguishable velocity 
components, whilst 12 are justified statistically. The present 
analysis shows the applicability of the comprehensive fitting 
method in a case of such a complex structure. However, it 
also shows that in order to achieve a competitive result on 
An/fi, the H2 absorber/spectrum selected for analysis must 
obey several conditions among which are no occurrence of 
strong additional absorbers and a sufficiently high signal-to- 
noise ratio. 
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